贝 叶 斯 结构 方程 模型 及 其 研究 现状 
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摘 要 在 心理 学 研究 中 结构 方程 模型 (Structural Equation Modeling，SEMD) 被 广泛 


速 ， 但 其 在 国内 心理 学 领域 的 应 


果 效 应 ， 其 估计 方法 有 频率 学 方法 (如 ， 极 大 似 然 估计 ) 和 贝 叶 斯 


及 其 在 结构 方程 建 模 中 易于 处 理 小 样本 、 


用 不 足 。 


常用 的 贝 叶 斯 结构 方程 模型 及 
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斯 估计 ; 极 大 似 然 估计 


1 引言 


在 心理 学 、 行为 学 和 其 他 社会 科学 


缺失 数据 及 复杂 模 


BR; WRZ” 


» 510006 


方法 两 类 。 近 年 来 


型 等 方面 的 优势 ， 贝 叶 斯 


1 叶 斯 结构 方程 模型 的 方法 基 


x 点 用 现状 ， =I 上 


究 领 域 , 研究 者 经 常 


LVs)， 即 不 能 直接 观测 的 变量 ,来 研 ? 
一 个 潜 变 量 往往 对 应 着 
作 其 对 应 外 显 变量 


特质 (Trait) 或 因子 (Factor)。 


特定 潜 变 量 的 反映 指标 。2 


析 洪 变量 间 关系 最 强 有 力 的 现代 统计 方法 , 是 心理 学 、 教 


方法 之 一 ( 修 杰 泰 ， 温 忠 鹿 


的 一 个 重要 原因 在 于 其 灵活 性 , 包括 
在 内 的 许多 模型 都 能 以 SEM 


传统 的 SEM 包括 测量 


EN 


(Manifest Variables)。 潜 变量 可 以 看 


介绍 新 的 研究 工具 。 


等 抽象 


用 于 检验 潜 变 量 间 的 


由 于 贝 叶 斯 统计 的 流行 


结构 方程 模型 发 展 迅 


础 和 优良 特性 ， 及 几 类 


要 通过 法 变量 (Latent Variables, 


的 构 念 (Construct)、 


F 个 相关 的 、 可 直接 观测 的 外 显 变量 


的 抽象 和 概括 ， 而 外 显 变量 则 可 视 为 


寺 构 方程 模型 (Structural Equation Modeling, SEM) 是 公认 的 用 于 分 


育 学 和 行为 


变量 关联 起 来 , 通常 是 一 个 验证 性 因 
于 检验 潜 变 量 的 因子 结构 。 其 
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定义 如 下 : 


量 模 型 将 外 显 变 量 与 对 应 的 潜 


学 等 领域 最 常用 的 分 析 
， 成 子 娟 ，2004; Lee & Song, 2012; Emi, 
验证 性 因子 分 析 模型 、 
的 形式 表征 
量 模型 和 结构 模型 两 个 部 分 ， 其 中 测量 


2014)。SEM 流行 


FP 介 模型 和 潜 变 量 增长 曲线 模型 


子 分 析 (Confirmatory Factor Analysis, CFA) 模 型 , 主要 用 


(1) 


上 的 观测 值 ，p (pX 1) 为 截 距 项 ，A (pXq) 


金 项 目 (18YJA190013); 中 山大 学 2018 


为 因子 载荷 矩阵 ， 用 于 反映 外 显 变量 六 和 潜 变 量 w (gq X 了 1) 之 间 的 关系 ，si;(pX1) 是 外 显 变 量 
的 测量 误差 项 ， 服 从 于 NI0, 亚 的 分 布 。 模 型 包含 以 下 假设 : 测量 误差 和 潜 变 量 间 不 存在 相 
关 ; 测量 误差 的 方差 协 方差 矩阵 于 ,为 对 角 和 矩 了 泗 ， 即 测量 误差 间 不 存在 相关 。 传 统 的 测量 模 
型 还 通常 假设 潜 变 量 的 数目 以 及 潜 变 量 与 外 显 变 量 间 的 关系 已 知 , 且 每 个 外 显 变量 只 负载 在 
一 个 潜 变 量 上 ， 不 存在 交叉 载荷 。 

结构 模型 主要 用 于 分 析 潜 变量 之 间 的 “因果 效应 ” 其 定义 如 下 : Sois ,6 )， 其 


中 i (qi X DESNE (Exogenous LVs), (qX 1) 2 A EVE (Endogenous LVs)， 


ni = Mn, + TE; + 6; = Aw; + 6; (2) 
其 中 AM, T), Mq X qi) TC X qo) ERE ARBOR, MRT Py AEP aS eZ Td A, 
而 工 则 反映 了 外 源 潜 变 量 对 内 生 潜 变量 的 效应 。& 服 从 于 MODIR, ôi (G1 X DAR 
项 ， 服 从 于 NOPR 
结构 方程 模型 的 分 析 步 又 通常 包括 模型 设 定 与 识别 、 模 型 拟 合 评 估 、 模 型 修正 和 参数 
估计 。 与 传统 的 路 径 分 析 方 法 相 比 ， 结 构 方程 模型 考虑 了 外 显 变量 的 测量 误差 ， 对 潜 变 量 间 
关系 的 估计 更 加 准确 ( 李 锡 钦 ，2011)。 


2 贝 叶 斯 结构 方程 模型 


结构 方程 模型 的 主要 估计 方法 包括 频率 学 派 方法 (如 ,， 极 大 似 然 估计 ) 和 贝 叶 斯 方法 两 类 。 
尽管 目前 对 频率 学 派 方法 的 应 用 更 多 , 但 近年 来 , 由 于 贝 叶 斯 方法 的 流行 及 其 在 统计 建 模 中 
= 的 诸多 优势 , 关于 贝 叶 斯 结构 方程 模型 的 方法 类 和 应 用 类 研究 数量 稳步 增长 。 尤 其 是 自 2012 


© 年 以 来 ， 贝 叶 斯 结构 方程 模型 的 应 用 研究 数量 大 幅 增 加 (Van de Schoot et al., 2017)。 


贝 叶 斯 方法 和 频率 学 派 方法 本 质 的 区 别 是 : 频率 学 派 将 未 知 参数 看 作 常 数 , 根据 样本 参 
数 估计 总 体 参 数 ; 而 贝 叶 斯 方法 则 将 未 知 参数 视 为 随机 变量 , 分 析 的 目的 是 得 到 未 知 参数 的 
后 验 分 布 ( 王 重 成 ， 邓 倩 文 ， 毕 向 阳 ，2017)。 在 用 贝 叶 斯 方法 分 析 SEM 时 ， 研 究 者 可 以 根 
据 理 论 或 以 往 研究 结果 确定 未 知 参 数 或 潜 变 量 的 先 验 分 布 , 如 果 没 有 准确 的 先 验 信息 也 可 以 
提供 无 信息 先 验 分 布 (如 ， 均 匀 分 布 ) 或 模糊 信息 先 验 分 布 (如 ， 方 差 极 大 的 正 态 分 布 ); 根据 
贝 叶 斯 公式 , 结合 先 验 分 布 和 数据 似 然 函数 可 以 得 到 未 知 参数 和 潜 变 量 的 后 验 分 布 ; 再 采用 
马尔 科 夫 链 蒙 特 卡 罗 (Markov Chain Monte Carlo，MCMC) 算 法 (如 ，Gibbs 抽样 法 和 
Metropolis-Hastings 算法 等 ) 从 后 验 分 布 中 从 代 地 抽取 大 量 样本 ; 通过 抽取 的 样本 估计 后 验 分 
布 的 均值 、 可 信 区 间 (Credible IntervaD 及 其 它 统计 量 ， 进 而 进行 统计 推断 ( 李 锡 钦 ，2011)。 


N 


采用 贝 叶 斯 方法 分 析 SEM 的 具体 步 又 包括 : 
(1) 设 定 模型 并 为 未 知 参 数 提供 先 验 信息 : S kale p, hale, q1， 对 于 SEM 中 不 同 


的 未 知 参数 提供 如 下 所 示 的 共 轿 先 验 分 布 ( 李 锡 饮 ，2011): 


-1 
H~N(Uo, Hpo), Ax~N(Aor, Hox), Aon ~N (Moon: Hown), © ~Wishart(Ro,Po) (3) 


其 中 A 是 载荷 矩阵 A 的 第 下行, Aor 是 路 径 系 数 和 矩阵 Ao 的 第 严 行 , po, Aow Aowns Po 

和 正定 矩阵 Ro, Ho, Hog Hoor 是 根据 理论 或 以 往 研究 结果 给 定 的 超 参数 值 , 反映 先 验 
信息 及 研究 者 对 先 验 信息 准确 性 的 把 握 ; 

(2) WE MCMC 算法 和 欠 代 次 数 ， 在 其 收敛 后 进行 模型 拟 合 评估 和 参数 估计 。 算 法 是 否 

达到 收敛 可 以 通过 自 相关 图 ， 踪 迹 图 和 潜在 尺度 缩减 因子 进行 评估 ( 详 见 王 孟 成 等 ， 
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y 2017); 
y (3) 模型 与 数据 的 整体 拟 合 程度 可 以 通过 后 验 预测 检验 (Posterior Predictive Checking) VP 
估 ; 


(4) 为 避免 先 验 信息 主观 性 的 影响 ， 研 究 者 可 以 通过 敏感 性 分 析 (Sensitivity Analysis; 
Greenland, 2001) 检 验 不 同 先 验 信息 下 估计 结果 是 否 稳定 ， 增 强 结果 的 可 靠 性 。 

与 传统 方法 相 比 ， 贝 叶 斯 结构 方程 模型 有 着 诸多 优势 ， 例 如 : (1) 相 比 于 频率 学 派 的 方 
法 ， 基 于 抽样 的 贝 叶 斯 方法 较 少 地 依赖 大 样本 渐 近 理论 ， 因 此 在 小 样本 中 依旧 表现 优良 
(Muthén & Asparouhov, 2012); (2) 贝 叶 斯 结构 方程 模型 的 分 析 基 于 原始 观测 值 ， 这 相 比 于 传 
统 方法 关注 的 协 方差 矩阵 更 易于 处 理 。 因 此 更 易于 处 理 复 杂 的 模型 和 数据 情况 , 如 存在 缺失 
一 值 的 数据 、 潜 变量 间 存 在 非 线性 关系 的 情况 等 ,而 传统 方法 在 这 种 情况 下 容易 遇 到 模型 识别 
© 问题 ( 李 锡 钦 ，2011); (3) 在 模型 拟 合 评估 、 模 型 比较 和 参数 估计 方面 ， 贝 叶 斯 方法 能 够 提供 

更 有 效 的 统计 量 (Pan, Ip, & Dubé, 2017); (4) 贝 叶 斯 方法 能 够 灵活 地 在 模型 估计 中 纳入 先 验 信 
息 ， 如 预 实 验 和 前 人 研究 结果 ， 而 有 效 的 先 验 信息 可 以 使 未 知 参 数 估计 更 加 准确 (Yuan & 
MacKinnon, 2009; Zhang, Hamagami, Wang, & Nesselroade, 2007). 


尽管 采用 贝 叶 斯 结构 方程 模型 有 着 诸多 优点 , 能 够 更 好 地 满足 应 用 研究 者 在 实证 研究 中 


的 需求 ， 如 处 理 复杂 模型 、 小 样本 问题 等 情况 , 但 其 在 国内 心理 学 领域 的 应 用 不 足 。 本 文 将 
介绍 几 类 常用 的 贝 叶 斯 结构 方程 模型 及 其 应 用 研究 进展 , 包括 验证 性 因子 分 析 、 结 构 方 程 模 
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3 分 类 


3.1 经 典 贝 叶 斯 结构 方程 模型 
3.1.1 贝 叶 斯 验证 性 因子 分 析 


验证 性 因子 分 析 常 被 用 于 根据 理论 假设 去 验证 外 显 变量 和 潜 变 量 间 的 关系 。 在 传统 的 

CFA 模型 中 ， 局 部 独立 性 (Local Independence) 假 设 要求 给 定 了 潜 变 量 的 值 后 ， 外 显 变量 之 间 

不 再 存在 相关 ， 即 测量 误差 项 的 方差 协 方差 矩阵 P, 中 非 对 角 线 元 素 的 值 均 被 限制 为 0。 此 

外 ， 传 统 的 CFA 模型 一 般 不 允许 存在 交叉 载荷 。 但 是 在 实际 研究 中 ， 这 种 严格 的 约束 条 件 

容易 导致 模型 拟 合 不 好 甚至 被 拒绝 (Muthén & Asparouhov, 2012)。 一 些 研究 者 指出 ， 传 统 方 

= 法 对 模型 施加 的 限制 是 过 于 严格 的 、 甚 至 是 不 必要 的 , 这 种 限制 在 大 样本 情况 下 很 容易 拒绝 

5 实际 上 和 数据 拟 合 良好 的 模型 (Lu,，Chow, & Loken, 2016; Marsh et al., 2009; Muthén & 

Asparouhov, 2012)。 且 有 研究 者 发 现 ， 在 实际 数据 分 析 中 对 模型 添加 过 多 的 约束 条 件 还 会 导 

致 未 知 参 数 估计 的 准确 
Thompson, 2014). 

在 传统 方法 中 ， 为 了 解决 这 种 限制 带 来 的 问题 ， 研 究 者 往往 会 结合 理论 和 修正 指数 
(Modification Index; Sérbom, 1989) 的 建议 ， 在 模型 中 增加 交叉 载荷 或 测量 误差 间 的 相关 。 但 
这 种 基于 修正 指数 的 方法 仍然 存在 一 些 局 限 ， 例 如 : (1) 由 于 需要 逐个 参数 进行 修正 ， 当 
要 修正 的 参数 较 多 时 ， 修 正 的 过 程 耗 时 、 繁 瑞 ，(2) 容 易 导 致 模型 的 过 拟 合 ， 削 弱 其 泛 化 


能 力 (Maccallum, Roznowski, & Necowitz, 1992); (3) 难 以 找到 全 局 最 优 的 模型 (Chou & Bentler, 


= 


生 降 低 (Asparouhov & Muthén, 2009; Hsu, Troncoso Skidmore, Li, & 


ad 1990); (4) 容 易 导 致 一 类 错误 率 增 大 (Draper 1995) 等 。 

Muthén 和 Asparouhov (2012) 创 新 地 提出 了 一 种 结合 了 探索 和 验证 方法 的 贝 叶 斯 验证 性 
因子 分 析 模 型 ， 放 宽 了 模型 对 于 测量 误差 相关 或 交叉 载荷 的 限制 。 在 传统 方法 中 ， 交 叉 载 荷 
和 测量 误差 间 的 相关 被 严格 限制 为 0， 这 既是 基于 模型 简洁 性 的 考虑 ， 也 是 因为 如 果 自 由 估 
计 这 些 参数 容易 导致 模型 不 可 识别 。 但 是 Muthén 和 Asparouhov (2012) 提 出 的 方法 在 保证 模 
型 可 识别 的 同时 ， 可 以 通过 对 交叉 载荷 提供 一 个 均值 为 0、 方差 极 小 的 正 态 先 验 分 布 ， 或 对 
误差 项 矩阵 提供 合适 的 道 Wishart 分 布 来 放宽 对 其 的 限制 。 模 拟 研究 显示 ， 该 方法 在 放宽 对 
于 交叉 载荷 或 测量 误差 相关 的 限制 时 , 得 到 的 显著 的 交叉 载荷 或 测量 误差 相关 的 数目 比 修正 
指数 方法 得 到 的 更 少 , 且 模 型 拟 合 在 一 次 分 析 中 就 可 以 得 到 满意 的 结果 , 而 传统 方法 通常 需 


要 进行 多 次 修正 。 


= 


但 是 Muthén 和 Asparouhov (2012) 提 出 的 这 种 方法 在 放宽 模型 限制 的 同时 ， 也 容易 导致 


较 多 非 零 交叉 载荷 或 测量 误差 相关 的 产生 (Lu, Chow, & Loken, 2016)， 使 得 因子 载荷 矩阵 或 
误差 项 矩阵 过 于 复杂 。 因 而 模型 容易 出 现 过 拟 合 的 情况 , 对 研究 结果 的 解释 和 重复 造成 困难 。 

Lu 等 人 (2016) 指 出 Muthén 和 Asparouhov (2012) 的 方法 本 质 上 是 将 贝 叶 斯 Ridge 正则 化 
(Regularization) 方 法 应 用 于 CFA 模型 中 。 针对 其 存在 的 上 述 问 题 , Lu 5 A(2016) 51 Af A 
种 贝 叶 斯 正则 化 方法 : 通过 对 载荷 矩阵 提供 spike-and-slab 先 验 分 布 , 保留 重要 的 交叉 载 答 ， 
将 其 它 微弱 的 交叉 载荷 压缩 到 零 。 这 种 方法 避免 了 Ridge 正则 化 方法 可 能 导致 的 模型 过 拟 合 ， 
及 其 对 重要 交叉 载荷 的 过 度 压缩 等 问题 。 

Pan 等 人 (2017) 则 针对 误差 项 的 方差 协 方差 矩阵 ， 将 协 方 差 Lasso (Least absolute 


shrinkage and selection operator) 正则 化 方法 引入 CEA 模型 。 通 过 估计 稀 疏 化 的 误差 协 方差 算 
阵 ， 在 放宽 对 测量 误差 相关 限制 的 同时 ， 将 微弱 的 、 不 重要 的 测量 误差 相关 向 零 压缩 ， 避 免 
因为 测量 误差 相关 过 多 而 导致 的 模型 过 拟 合 或 误差 项 矩阵 不 正定 等 问题 。 其 实证 研究 发 现在 
允许 “少量 ”测量 误差 相关 的 情况 下 ，CFA 模型 的 简约 性 和 拟 合 程度 都 得 到 了 满足 。 

自 放宽 模型 限制 的 思想 被 提出 以 来 ， 由 于 其 诸多 优点 ， 采 用 贝 叶 斯 CFA 进行 数据 分 析 
的 应 用 研究 越 来 越 多 。 基 于 贝 叶 斯 方法 在 处 理 复杂 模型 时 的 优良 特性 ，Golay，Reverte， 
Rossier, Favez 和 Lecerf (2013) 重 新 分 析 了 韦 氏 智力 量 表 的 四 因子 结构 , 分 别 检验 了 二 阶 因子 
模型 和 双 因 子 (Bifactor) 模 型 , 结果 显示 贝 叶 斯 方法 在 模型 识别 和 估计 上 都 比 极 大 似 然 估计 表 
现 得 更 好 ;而 Falkenstrom 等 人 (2015) 在 检验 病人 版 工作 智力 量 表 的 结构 效 度 时 发 现 ， 极 大 
似 然 估计 显示 模型 拟 合 较 差 , 但 采用 贝 叶 斯 方法 放宽 对 测量 误差 相关 的 限制 后 , 模型 和 数据 
拟 合 很 好 ， 此 外 ， 由 于 小 样本 中 贝 叶 斯 方法 对 于 参数 的 估计 更 加 准确 (Muthén & Asparouhov, 


2012)，Crenshaw, Christensen, Baucom, Epstein 和 Baucom (2016) 在 临床 的 小 样本 研究 中 使 用 
贝 叶 斯 CFA 方法 修订 了 沟通 模式 量 表 (52 名 被 试 ， 模 型 有 18 个 未 知 参数 ， 包 括 9 个 因子 载 
荷 和 9 个 测量 误差 方差 ); 而 由 于 与 传统 的 极 大 似 然 估计 方法 相 比 ， 贝 叶 斯 方法 在 小 样本 的 
情况 下 对 因子 分 的 估计 更 加 准确 (Muthén & Asparouhov, 2012), Alessandri 和 De Pascalis 
(2017) 在 脑 电 实验 研究 中 采用 贝 叶 斯 CFA 估计 51 名 被 试 “生活 导向 ”因子 的 因子 分 ， 再 用 
于 后 续 因 子 间 关系 的 分 析 中 。 


a 


3.1.2 贝 叶 斯 结构 方程 模型 


在 测量 模型 的 基础 上 ,结构 模型 被 用 于 检验 不 受 测量 误差 影响 的 潜 变 量 间 的 “因果 效应 ”。 
正如 本 文 第 二 部 分 所 述 , 在 采用 贝 叶 斯 方法 估计 结构 方程 模型 时 , 可 以 通过 对 未 知 参数 和 洪 
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变量 提供 合适 的 


先 验 分 布 ， 结 合 数 据 得 到 其 后 验 分 布 ， 进 而 进行 统计 推断 。 


随 着 对 测量 


在 建立 测量 模型 


验 分 布 。 和 贝 叶 斯 验证 性 因子 分 析 相 似 的 是 ，Muthén 和 Asparouhov (2012) 指 出 


模型 限制 的 放宽 , 对 结构 模型 中 未 知 参数 的 估计 也 更 加 准确 。 Pan 等 人 (2017) 


时 采用 贝 叶 斯 Lasso 方法 对 误差 协 方差 矩阵 进行 估计 ， 结 果 发 现 : 与 传统 方 
法 相 比 ， 贝 叶 斯 Lasso 方法 下 对 结构 模型 中 路 径 系数 的 估计 偏差 也 更 小 。 
此 外 ， 在 对 结构 模型 的 估计 中 ， 使 用 贝 叶 斯 方法 可 以 对 所 有 的 路 径 系数 提 


供 一 定 的 先 
可 以 对 原先 


被 限定 为 0 的 路 径 系 数 提供 均值 为 0， 方差 极 小 的 正 态 先 验 分 布 来 放宽 对 其 的 限制 ， 在 一 次 


估计 中 可 以 同时 
为 了 方便 


构 方程 模型 为 例 ， 


研究 重新 分 析 了 
所 建立 的 模型 如 
过 影响 其 十 年 


实现 对 模型 探索 和 验证 。 


究 者 应 用 这 种 方法 , 我们 将 以 Muthén 和 Asparouhov (2012) 建 立 的 贝 叶 斯 结 


详细 介绍 在 测量 模型 和 结构 模型 中 同时 使 用 信息 先 验 的 具体 实现 步 又 。 访 


Kaplan (2009) 的 追踪 研究 ,研究 收集 了 6677 名 公立 学 校 十 年 级 学 生 的 数据 ， 


图 1 所 示 。 模 型 包含 以 下 假设 : 学 生 六 年 级 的 科学 成 绩 和 家 庭 社 
级 的 科学 成 绩 进而 影响 学 生 的 科学 成 就 ; 而 教师 的 科学 认证 水 平 通 


生 的 参与 程度 , 3 


响 其 科学 成 就 。 


六 年 级 科学 


影响 学 生 对 课程 挑战 性 的 感知 程度 , 最 终 通过 学 生 十 年 级 的 


成 绩 十 年 级 科学 成 绩 科学 成 就 


家 庭 社会 经 济 地 位 


会 经 济 地 位 
过 影响 学 


科学 成 绩 影 


教师 科学 


学 认证 水 平 


观察 观察 观察 观察 观察 观察 
指标 1 指标 2 指标 3 指标 4 指标 5 指标 6 


研究 者 采用 传统 的 极 大 似 然 估计 发 现 ， 


图 1 科学 成 就 的 结构 方程 模型 图 (改编 自 Kaplan(2009)) 


CFI=0.844)， 修 正 指数 显示 模型 中 存在 大 量 需 要 被 修正 的 参数 。 而 无 信息 先 验 下 


的 估计 结果 与 之 


相似 ， 后 验 预测 检验 区 间 不 包括 零 ， 模 型 拟 合 较 差 (959%CI[1644, 
基于 前 文中 放宽 模型 限制 的 思想 ， 在 对 数据 进行 标准 化 处 理 后 ，Muthen 和 


(2012) 针 对 结构 模型 中 原先 被 固定 为 0 的 11 个 参数 提供 了 均值 0， 方 差 为 0.01 


模型 拟 合 较 差 , 假设 模型 被 拒绝 (RMSEA=0.081, 


贝 叶 斯 方法 
1720])。 
Asparouhov 


的 正 态 分 布 


先 验 ， 即 允许 这 些 参数 有 95% 的 变异 落 在 (-0.2,，0.2) 之 间 。 由 于 已 经 对 数据 进行 了 标准 化 处 
理 ， 上述 设 定 的 先 验 分 布 理 论 上 能 较 好 地 覆盖 那些 可 能 被 忽略 的 效应 的 范围 。 研究 者 针对 六 
年 级 科学 成 绩 、 家 庭 社会 经 济 地 位 和 教师 科学 认证 水 平 三 个 变量 指向 六 个 观察 指标 的 18 个 
路 径 系数 也 提供 了 同样 的 先 验 分 布 , 在 传统 方法 中 , 这 些 路 径 系数 如 果 被 自由 估计 将 意味 着 
六 个 观察 指标 存在 项 目 功 能 差异 。 但 是 在 这 里 ， 研 究 者 仅仅 通过 提供 均值 为 0， 方 差 极 小 的 
先 验 分 布 来 放宽 对 其 严格 的 限制 。 

此 外 ， 研 究 者 还 放宽 了 对 于 六 个 观察 指标 的 交叉 载荷 及 测量 误差 相关 的 限制 。 其 中 ， 对 
交叉 载荷 同样 提供 了 均值 为 0, 方差 为 0.01 的 正 态 先 验 分 布 ， 而 对 误差 项 矩阵 提供 IWA, 30) 
fit Wishart 分 布 。 模 型 对 应 的 Mplus 代码 详 见 Muthén 和 Asparouhov (2012) 的 补充 材料 。 

结果 发 现 与 未 放宽 模型 限制 的 模型 相 比 ， 模 型 拟 合 得 到 了 提升 ， 后 验 预 测 检验 区 间 包 
括 0 (95% CI [-24, 44]), 后 验 预测 疡 值 为 0.276。 同 时 还 额外 发 现 了 6 个 显著 的 结构 模型 参数 ， 
及 11 对 显著 的 测量 误差 相关 ， 但 其 估计 值 均 较 小 !。Muthén 和 Asparouhov (2012) 在 文章 中 
指出 ， 如 果 使 用 传统 的 极 大 似 然 估计 ， 上 述 放宽 限制 后 的 模型 将 是 不 可 识别 的 , 但 在 贝 叶 斯 
框架 下 可 以 避免 这 种 情况 。 同 时 他 们 还 指出 在 放宽 模型 的 限制 后 ， 模 型 的 收敛 率 降低 ， 需 要 
更 多 的 迭代 次 数 和 更 长 的 运算 时 间 才 能 得 到 模型 估计 结果 。 关 于 上 述 内 容 的 讨论 , 感 兴趣 的 
读者 可 参考 Muthén 和 Asparouhov (2012), A ICA FLEA. 

此 外 ， 在 实证 研究 中 ， 越 来 越 多 的 应 用 研究 者 也 开始 采用 这 种 方法 。Scherer, Siddiq 和 
Teo (2015) 基 于 所 研究 的 潜 变 量 间 存 在 概念 重合 的 情况 ， 认 为 建 模 时 必须 要 考虑 交叉 载 答 ， 
但 传统 方法 中 加 入 过 多 的 交叉 载荷 会 出 现 模型 识别 或 估计 问题 基于 贝 叶 斯 结构 方程 模型 的 
PE 思想 ，Scherer  A(2015) AGFA T THII AA EA ERR R, JEE EEN 

模型 以 检验 教师 对 信息 技术 的 利用 与 教师 自身 特征 的 关系 。Salarzadeh, Moghavvemi, Wan,， 
Babashamsi 和 Arashi (2017) 同 样 采用 该 方法 研究 了 影响 学 生 使 用 电子 学 习 平台 意愿 的 因素 ， 
并 对 比 极 大 似 然 估计 和 贝 叶 斯 估计 的 结果 , 发 现 贝 叶 斯 估计 下 模型 估计 结果 更 准确 , 误差 均 
方 根 和 绝对 平均 误差 更 小 。 


co 


上 建立 结构 


3.1.3 贝 叶 斯 中 介 模 型 


I 


中 介 分 析 在 心理 学 领域 中 扮演 着 非常 重要 的 角色 , 它 可 以 用 于 解释 自 变 量 对 因 变 量 的 作 
用 机 制 ， 有 助 于 已 有 理论 的 验证 和 新 的 理论 的 构建 ( 罗 胜 强 ， 姜 录 ，2014)。 在 检验 中 介 效 应 


lin 


! 研究 者 将 显著 的 测量 误差 相关 仅 视 为 元 余 参数 (nujisance parameters). 


时 ， 关 键 是 检验 间接 效应 的 显著 性 ， 即 自 变量 是 否 会 通过 中 介 变 量 “ 显 著 地 ”影响 因 变 量 。 


在 Yuan 和 MacKinnon (2009) 的 贝 叶 斯 中 介 建 模 方 法 出 现 之 前 ， 心 理学 研究 中 大 多 数 中 介 效 


应 分 析 都 在 频率 学 派 的 框架 下 进行 (Nuijten, Wetzels, Matzke, Dolan, & Wagenmakers, 2015). 


目前 频率 学 派 中 检验 间接 效应 显著 性 的 常用 方法 是 Sobel 法 和 Bootstrap 法 。 其 中 ，Sobel 法 
即 直 接 检 验 直接 效应 系数 乘积 的 显著 性 (Hu: ab=0; a, b 即 两 个 直接 效应 系数 )， 但 是 该 方法 
需要 假设 系数 的 乘积 服从 正 态 分 布 , 否则 估计 将 产生 偏差 。 而 这 个 假设 在 实际 数据 分 析 中 往 
往 很 难 被 满足 ,但 Bootstrap 法 通过 构建 区 间 估 计 可 以 避免 这 一 问题 ( 温 忠 腾 , 叶 宝 娟 ,2014)。 

在 贝 叶 斯 框架 下 的 中 介 分 析 则 是 通过 MCMC 法 进行 ， 这 种 方法 基于 从 后 验 分 布 中 抽取 
的 样本 进行 参数 估计 。 在 获取 了 直接 效应 系数 的 后 验 分 布 后 ， 可 以 很 容易 地 对 参数 的 各 种 
函数 形式 进行 估计 ， 如 两 个 直接 效应 系数 的 乘积 ， 而 且 该 方法 易于 构造 间接 效应 的 可 信 区 
间 (Yuan & MacKinnon, 2009)。 因 此 贝 叶 斯 方法 易于 处 理 复杂 的 中 介 模 型 ， 如 贝 叶 斯 序列 中 


x 


介 模 型 (Tofighi & Mackinnon, 2016)， 有 调节 的 中 介 模 型 (Wang & Preacher, 2015) 等 。 而 传统 


方法 在 处 理 复杂 模型 时 则 常常 会 遇 到 模型 识别 和 估计 问题 (Kenny, Korchmaros, & Bolger, 


2003)。 

此 外 ， 贝 叶 斯 中 介 模 型 不 需要 假设 系数 乘积 服从 正 态 分 布 。 以 往 的 研究 发 现 ， 除 无 信 
息 先 验 条 件 下 贝 叶 斯 方法 和 传统 方法 的 估计 结果 比较 相似 外 ， 在 信息 先 验 和 模糊 信息 先 验 
条 件 下 ， 贝 叶 斯 方法 的 检验 力 都 更 高 ， 对 于 参数 的 估计 也 更 加 准确 (MacKinnon, Lockwood, 


& Williams, 2004; Tofighi & MacKinnon, 2011; Tofighi & Mackinnon, 2016)， 在 小 样本 中 表现 


也 更 好 (Miotevié, MacKinnon, & Levy, 2017). 

由 于 中 介 分 析 的 重要 性 及 贝 叶 斯 中 介 模 型 的 优良 特性 ， 采 用 贝 叶 斯 方法 检验 中 介 效 应 
的 研究 也 越 来 越 多 。 鉴 于 贝 叶 斯 分 析 不 依赖 于 参数 的 正 态 分 布 假设 , Shuck, Zigarmi 和 Owen 
(2015) 通 过 贝 叶 斯 多 重 中 介 模 型 检验 了 工作 投入 、 员 工 敬业 度 、 工 作 热 情 在 基本 心理 需求 对 
工作 意图 的 影响 中 的 中 介 效 应 ， 而 由 于 研究 的 样本 量 较 少 ，Zeman，Dallaire，Folk 和 Thrash 
(2017) 同 样 采 用 该 方法 检验 了 在 监禁 风险 经 历 和 环境 风险 对 儿童 精神 障碍 的 影响 中 儿童 情 
绪 管 理 的 中 介 效 应 ， 此 外 ， 结 合 贝 叶 斯 中 介 建 模 和 贝 叶 斯 CFA 放宽 测量 模型 限制 的 方法 ， 


Jacobson, Lord 和 Newman (2017) 验 证 了 焦虑 通过 情绪 智力 中 介 影 响 抑 郁 症状 。 


3.2 贝 叶 斯 潜 变 量 增 长 曲线 模型 


纵向 研究 也 称 追 踪 研 究 , 它 通 过 对 研究 对 象 进行 长 期 追踪 观察 .重复 测量 个 体 有 关 变 量 ， 


能 够 观察 较为 完整 的 发 展 过 程 ， 发 现 发 展 过 程 中 的 一 些 关 键 转 折 。 与 横断 研究 相 比 ， 纵 向 研 


究 设计 最 大 的 优点 是 可 以 合理 地 推导 变量 之 间 的 因果 关系 , 是 心理 学 研究 中 一 种 重要 的 研究 
方法 ( 刘 红 云 ， 重 庆 成 ，2003)。 
目前 在 分 析 追 踪 数 据 的 结构 方程 建 模 中 ， 潜 变量 增长 曲线 模型 (Latent Growth Curve 
Model, LGCM; Bollen & Curran, 2006) 应 用 非常 广泛 。 该 模型 不 仅 能 够 分 析 总 体 平 均 变化 趋 
势 ， 还 能 分 析 不 同 个 体 发 展 轨迹 变化 的 差异 及 其 如 何 受 到 预测 变量 的 影响 。 

LGCM 适用 于 在 几 个 固定 时 间 点 观测 得 到 的 纵向 研究 数据 ( 刘 红 云 ， 备 庆 成 ，2003)。 一 
个 简单 的 非 条 件 线性 潜 变 量 增 长 曲线 模型 的 定义 如 下 ( 王 济 川 ， 王 小 倩 ， 姜 宝 法 ，2011): 


yti = Noi + Atai + Eti (4) 
Noi = No + Goi (5) 
mi=mt oii (6) 


HP yrth i PAREA tT, go 和 辐 是 随机 系数 ， 描 述 了 数据 随时 间 
变化 的 特点 ， 也 被 称 为 发 展 因子 (Growth Factors)， 久 是 时 间 分 值 ，si; 是 变量 在 时 间 点 1 的 复 
合 误 差 项 ， 表 示 随 机 测量 误差 和 第 i 个 个 体 的 特定 时 间 效 应 。 在 (5) 式 和 (6) 式 中 ，7Jo 和 分 
别 反 映 变 量 的 平均 初始 水 平和 变量 的 平均 变化 率 , 故而 yo 和 jy1; 分 别 被 称 为 潜 截 距 发 展 因 子 
AUPE RR ARAL» Co 和 61; 表示 残 差 项 ， 其 方差 能 反映 研究 对 象 的 个 体 间 变异 。 

此 外 , 研究 者 还 可 以 在 模型 中 纳入 一 些 预测 变量 来 预测 上 述 两 个 潜 发 展 因 子 。 而 如 果 变 
量 的 发 展 轨 迹 为 非 线 性 , 也 可 以 通过 引入 时 间 分 值 的 三 阶 或 高 阶 函 数 , 或 通过 设 定 自由 时 间 
分 值 等 方式 来 检验 非 线 性 的 发 展 轨 迹 (Meredith & Tisak, 1990)。 

在 估计 潜 变 量 增长 曲线 模型 时 ， 贝 叶 斯 方法 相 比 于 传统 估计 方法 上 共有 许多 优势 。 首 先 ， 
追踪 数据 情况 较为 复杂 ， 可 能 会 出 现 测量 时 间 间 隔 不 相等 、 发 展 轨迹 非 线 性 、 数 据 不 满足 正 
态 分 布 等 情况 , 因此 有 时 需要 进行 复杂 的 数据 转换 ， 而 贝 叶 斯 估计 依赖 于 原始 观测 值 而 非 方 
差 协 方差 矩阵 ， 更 易于 处 理 这 种 数据 转换 的 情况 (Gelman，Carlin，Sterm，Dunson，Vehtari， 人 


Rubin, 2014); 其 次 ， 传 统 方法 在 处 理 复杂 的 LGCM 时 容易 出 现 估计 问题 甚至 是 模型 识别 问 


题 ， 而 贝 叶 斯 方法 能 够 更 有 效 地 估计 复杂 模型 (Zhang, Hamagami, Wang, & Nesselroade, 2007); 
第 三 ， 追 踪 研 究 中 普遍 存在 缺失 数据 ( 叶 素 静 ， 唐 文清 ， 张 敏 强 ， 曹 魏 聪 ，2014)， 而 贝 叶 斯 
方法 在 处 理 缺 失 数据 时 表现 更 好 ( 李 锡 钦 ，2011D; 此 外 ， 贝 叶 斯 方法 在 小 样本 情况 下 表现 更 
好 (Gelman, Carlin, Stern, & Rubin, 2003), Zhang 等 人 (2007) 发 现 使 用 贝 叶 斯 方法 甚至 可 以 估 
计 只 有 20 个 被 试 在 4 个 不 同时 间 点 测量 的 潜 变 量 增长 曲线 模型 ， 这 使 得 贝 叶 斯 方法 具有 极 
高 的 应 用 价值 。 

在 实证 研究 中 ，Winans-Mitrik 等 人 (2014) 指 出 比 起 依赖 假设 检验 的 传统 频率 学 派 方法 ， 


贝 叶 斯 方法 能 够 提供 更 加 稳定 的 推断 ， 更 符合 临床 评估 的 要 求 ， 因 而 采用 贝 叶 斯 LGCM 验 
证 失语 症 干预 项 目 对 于 语言 理解 的 促进 作用 , 在 对 被 试 的 失语 症 严 重 程度 进行 了 四 次 测量 后 ， 
发 现 该 干预 项 目 对 于 失语 症 有 可 靠 且 持续 的 治疗 作用 。 此外， 由 于 贝 叶 斯 方法 在 小 样本 和 数 
据 非 正 态 情 况 下 表现 更 优良 (Gelman et al., 2014), Maier, Bohlmann 和 Palacios (2016) 采 用 贝 
叶 斯 LGCM 分 析 了 儿童 语言 发 展 中 跨 语 言 的 联系 ， 发 现 英语 与 西班牙 语 的 词汇 表达 能 力 的 
交互 项 影响 英语 词汇 表达 能 力 的 增长 。 我 们 将 以 这 个 研究 为 例 ， 详 细 介 绍 贝 叶 斯 LGCM 的 
建 模 步 又 。 

该 研究 的 被 试 为 洛杉矶 地 区 低 到 中 等 收入 家 庭 中 的 177 名 学 龄 前 儿童 ,研究 者 在 秋季 学 
期 初 、 四 个 月 后 、 六 个 月 后 三 个 时 间 点 对 其 进行 随访 。 在 每 次 随访 中 ， 通 过 语言 测试 及 编码 
员 进 行 行为 编码 的 方式 , 测量 儿童 英语 和 西班牙 语 的 词汇 表达 能 力 与 接受 能 力 , 并 调查 了 儿 


> 童 的 人 口 学 变量 、 教 师 与 班级 特征 、 班 级 与 家 庭 语言 环境 等 方面 的 多 个 协 变量 。 

© 研究 者 通过 贝 叶 斯 方法 处 理 缺 失 值 , 将 缺失 观测 值 作为 未 知 值 进行 估计 。 在 进行 建 模 分 

= 析 时 ， 首 先 为 英语 词汇 接受 能 力 、 英 语词 汇 表 达能 力 、 西 班 牙 语词 汇 接受 能 力 、 西 班 牙 语词 

N 汇 表达 能 力 四 个 结果 变量 分 别 建立 非 条 件 潜 变量 增长 曲线 模型 (图 2). 

CO 模型 为 每 个 变量 构造 了 两 个 发 展 因子 , 截 距 与 斜率 因子 的 协 方差 被 允许 自由 估计 。 每 个 

= 因子 都 有 对 应 三 个 时 间 点 的 测量 作为 它们 的 观察 指标 .对 于 斜率 因子 , 根据 测量 的 时 间 间 陋 ， 
按照 线性 关系 固定 了 三 次 时 间 点 的 负载 分 别 为 0、0.4、1。 而 三 个 时 间 点 中 观察 指标 的 残 关 


被 允许 自由 估计 。 截 距 和 斜率 发 展 因子 的 均值 反映 了 该 变量 的 总 体 平 均 初 始 值 及 变化 率 , 因 
子 的 方差 则 反映 了 个 体 间 差异 。 在 模型 估计 时, 研究 者 采用 了 Mplus 6.11 中 的 贝 叶 斯 估计 并 
提供 默认 的 无 信息 先 验 分 布 。 而 在 建立 了 非 条 件 模型 之 后 ,研究 者 进一步 将 协 变量 纳入 到 模 
型 中 以 预测 截 距 因子 和 斜率 因子 ， 建 立 条 件 潜 变 量 增长 曲线 模型 。 

结果 发 现 被 试 群体 的 英语 词汇 接受 能 力 、 英 语词 汇 表 达能 力 、 西 班 牙 语 词汇 接受 能 力 以 
及 西班牙 语词 汇 表 达能 力 对 应 的 四 个 非 条 件 模 型 拟 合 良 好 ，PPp 值 分 别 为 0.56( 后 验 预 测 检 
验 区 间 95%CI [-13.51, 10.46]), 0.26(95%CI [-8.38, 16.85]), 0.53(95%CI [-11.33，9.66])， 
0.48(95%CI [-12.67, 11.81])。 各 个 斜率 因子 估计 结果 均 大 于 0， 表明 四 种 词汇 能 力 均 随 着 时 
间 而 逐渐 增长 。 在 四 种 能 力 上 ， 不 同 个 体 在 初始 水 平和 变化 率 上 表现 出 显著 的 个 体 差 异 ( 即 
研 距 和 和 斜率 因子 的 方差 均 显 著 不 为 0)。 


加 入 协 变量 后 的 条 件 模 型 的 PPp 值 分 别 为 0.31(95%CI [-26.96, 52.53]), 0.21(95%CI 


[-23.84, 58.95]), 0.29(95%CI [-32.40, 51.26]), 0.34(95%CI [-28.35, 52.86]) ， 模 型 拟 合 结果 均 


为 良好 。 参 数 估计 结果 显示 ， 在 四 个 条 件 模型 中 ,年 龄 对 截 距 因 子 都 有 正 向 预测 作用 ， 表 明 
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年 龄 较 大 的 学 龄 前 儿童 倾向 于 拥有 较 高 的 初始 词 》 


量 。 家 庭 中 的 英语 接触 水 平 也 能 够 正 向 预 


测 儿 童 英语 词汇 能 力 ， 而 对 西班牙 语 的 初始 词汇 能 力 有 负 向 影响 。 另 外 ， 由 于 贝 叶 斯 方法 在 
估计 因子 得 分 上 的 优势 ， 研 究 者 在 贝 叶 斯 LGCM 下 得 到 了 每 个 被 试 在 截 距 因 子 和 和 斜率 因子 
上 的 值 ， 进一步 探究 了 英语 和 西班牙 语词 汇 能 力 间 的 相互 联系 ,但 由 于 本 文 篇 幅 有 限 ， 在 此 


BURR o 


英语 词汇 接受 
能 力 
截 距 因 子 


英语 词汇 英语 词汇 
接受 能 力 接受 能 力 


Tl T2 


英语 词汇 
接受 能 力 
T3 


图 2 非 条 件 潜 变 量 增长 曲线 模型 (以 英语 词汇 接受 能 力 为 例 ) 
上 述 实例 简单 介绍 了 贝 叶 斯 方法 在 潜 变 量 增 长 曲线 模型 中 的 应 用 。 目 前 , 在 大 部 分 采用 
贝 叶 斯 LGCM 的 研究 中 ， 研 究 者 仅仅 是 将 模型 估计 方法 改 为 贝 叶 斯 估计 ， 对 于 参数 估计 往 


往 提 供 软 件 默 认 的 无 信息 先 验 。 未 来 研究 可 以 考察 合适 的 先 验 分 布 对 于 潜 变 量 增长 曲线 模型 
中 参数 估计 的 影响 ， 进 一 步 发 挥 贝 叶 斯 方法 在 模型 拟 合 和 参数 估计 中 的 优势 。 


3.3 贝 叶 斯 多 组 结构 方程 模型 


在 实证 研究 中 , 研究 者 常常 会 根据 某 些 变 量 将 人 群 分 为 多 个 组 别 , 如 根据 性 别 划分 出 男 
女 两 组 。 在 这 种 情况 下 ,数据 往往 表现 为 组 数 较 少 、 
的 ， 即 多 组 数据 。 在 多 组 数据 的 结构 方程 建 模 中 , 看 


每 组 观测 值 较 多 且 组 内 观测 值 均 是 独立 


究 者 可 以 通过 分 组 建立 模型 来 研究 不 同 


组 间 模 型 的 相似 性 和 差异 性 ,多 组 结构 方程 模型 所 感 兴趣 的 是 检验 不 同 组 间 模 型 的 各 种 不 变 


性 假设 ( 李 锡 钦 ，2011)。 


不 变性 假设 的 检验 首先 


需 
假设 。 检 验 测量 不 变性 可 以 通过 一 系列 峰 套 步骤 进行 ， 按 顺序 包括 结构 不 变性 (各 组 具有 相 


要 在 测量 模型 中 进行 ， 即 测量 不 变性 (Measurement Invariance) 


同 的 因子 结构 )， 即 量 表 条 目 在 不 同 组 中 均 被 负载 到 相同 的 潜 变 量 上 ， 载 荷 不 变性 (条 目 载荷 


跨 组 相同 )， 截 距 不 变性 (条 目 截 距 跨 组 相同 ) 以 及 误差 方差 / 协 方差 不 变性 (条 目 误 差 方 差 跨 组 
相同 )(Vandenberg & Lance, 2000). 测量 不 变性 是 进行 潜 变 量 水 平 参 数 的 跨 组 比较 的 基础 ， 即 
对 结构 不 变性 的 检验 ， 结 构 不 变性 包括 因子 方差 / 协 方差 不 变性 和 因子 均 数 结构 不 变性 ( 王 济 
川 等 ，2011)。 当 模型 不 符合 跨 组 截 距 不 变性 时 ， 跨 组 因子 分 的 差异 通常 不 能 清楚 地 反映 潜 
变量 水 平 的 真实 差异 (Schmitt & Kuljanin, 2008)。 

在 检验 不 变性 时 ， 传 统 方法 首先 需要 确定 每 个 组 的 基线 模型 ， 然 后 构建 组 态 模型 
(Configural Model)， 即 整合 各 组 的 基线 模型 ， 在 组 态 模型 中 对 相应 参数 (如 ， 载 苟 、 截 距 ) 施 
加 跨 组 不 变 的 限制 ; 再 根据 模型 拟 合 变化 情况 判断 是 否 存在 违背 跨 组 不 变 的 参数 , 如 果 模 型 
拟 合 变化 显著 ， 则 需要 根据 修正 指数 的 建议 进行 参数 修正 ( 王 济 川 等 ，2011)。 当 违背 不 变性 
的 参数 占 被 限制 跨 组 不 变 的 参数 数目 比重 较 小 时 ， 这 种 修正 对 参数 估计 的 准确 性 影响 不 大 。 
但 是 当 组 数 较 多 时 , 违背 不 变性 的 参数 通常 较 多 , 这 种 情况 下 模型 修正 容易 导致 有 偏 的 参数 
估计 。 此 外 ， 传 统 方法 需要 根据 修正 指数 逐个 参数 进行 修正 ， 建 模 过 程 繁琐 ， 模 型 复杂 ， 而 
且 常 常 造成 模型 拟 合 较 差 和 参数 的 有 偏 估计 (Asparouhov & Muthén, 2014). 

而 贝 叶 斯 方法 可 以 很 好 地 处 理 组 数 较 多 的 情况 。 它 通过 对 参数 跨 组 的 差异 值 提供 合适 的 
先 验 分 布 (如 ， 均 值 为 0、 方 差 极 小 的 正 态 分 布 ) 来 放宽 对 参数 严格 跨 组 不 变 的 限制 ， 一 定 程 
度 上 人 允许 参数 跨 组 存在 较 小 的 差异 ， 从 而 避免 限制 过 于 严格 而 导致 的 模型 拟 合 过 差 的 问题 ， 
同时 也 可 以 提供 更 准确 的 参数 估计 (Asparouhov & Muthén, 2014)。 此 外 ， 贝 叶 斯 方法 在 处 理 
小 样本 的 多 组 数据 时 表现 也 更 好 (Kim et al., 2017). 

考虑 到 传统 方法 对 模型 施加 了 过 多 的 限制 , 容易 导致 模型 拟 合 较 差 。 Fong (2014) 结 合 贝 
叶 斯 多 组 建 模 和 贝 叶 斯 CFA 放宽 测量 模型 限制 的 方法 ， 检 验 了 男女 两 组 在 医院 焦虑 抑郁 
(Hospital Anxiety and Depression) 得 分 上 的 跨 组 不 变性 , 结果 发 现 女 性 焦虑 水 平 显 著 高 于 男性 。 


而 De Bondt 和 Van Petegem (2015) 在 该 方法 的 基础 上 ,基于 Asparouhov 和 Muthén (2014) 
的 建议 , 通过 对 参数 跨 组 的 差异 值 提供 先 验 分 布 来 放宽 对 于 参数 严格 跨 组 不 变 的 限制 , 检验 
过 度 激 动 问卷 的 性 别 跨 组 不 变性 。 发 现在 传统 方法 中 模型 拟 合 被 拒绝 的 多 组 模型 , 在 放宽 
对 参数 的 严格 限制 后 ， 拟 合 良 好 ， 满 足 了 测量 不 变性 。 我 们 将 以 该 研究 为 例 ， 详 细 介 绍 贝 叶 
斯 多 组 结构 方程 建 模 的 具体 实现 步骤 。 

De Bondt 和 Van Petegem (2015) 采 用 贝 叶 斯 多 组 结构 方程 模型 来 检验 过 度 激动 问卷 的 
别 跨 组 不 变性 。 研 究 收集 了 516 名 大 学 生 的 数据 ， 要求 被 试 在 网 上 填写 过 度 激动 问卷 ， 其 中 
问卷 包括 精神 运动 、 感 官 、 智 力 、 想 象 和 情感 5 个 维度 ， 共 50 题 。 

研究 者 首先 分 别 建立 了 男女 两 组 的 测量 模型 , 估计 方法 分 别 采 用 了 传统 的 频率 学 派 和 贝 


f= 


叶 斯 方法 。 在 贝 叶 斯 方法 中 ,放宽 了 对 于 测量 误差 相关 和 交叉 载荷 的 限制 : 其 中 对 交叉 载荷 
提供 均值 为 0, 方差 为 0.01 的 正 态 分 布 先 验 ; 而 对 观察 指标 的 误差 协 方差 矩阵 提供 IWA, 56) 
的 逆 Wishart 先 验 分 布 。 结 果 显 示 : 传统 方法 下 男女 两 组 相应 模型 的 拟 合 均 较 差 (CFI 值 均 小 
于 0.8), 而 贝 叶 斯 方法 下 男女 两 组 的 测量 模型 拟 合 良 好 , 模型 PPp 值 分 别 为 0.905 和 0.767, 
95% 区 间 包 括 零 ， 满 足 模型 跨 组 结构 不 变性 。 

在 此 基础 上 , 研究 者 建立 了 五 个 多 组 模型 ,比较 了 男女 性 在 过 度 激 动 问卷 5 个 维度 上 的 
跨 组 不 变性 。 作 为 参照 组 ,研究 者 将 男性 组 别 中 的 因子 均值 和 方差 则 分 别 固定 为 0 和 1， 而 
对 女性 组 别 中 的 因子 均值 和 方差 , 以 及 男女 组 间 的 协 方差 提供 无 信息 先 验 分 布 , 即 均值 为 0， 
方差 极 大 的 正 态 先 验 分 布 ， 即 进行 自由 估计 以 比较 组 间 因 子 均值 差异 。 此 外 ,研究 者 不 仅 放 
宽 了 对 于 测量 误差 相关 和 交叉 载荷 的 严格 限制 , 还 通过 对 因子 载荷 和 截 距 的 组 间 差 异 值 提供 
均值 为 0, 方差 为 0.01 的 正 态 分 布 先 验 ， 放宽 对 载荷 不 变性 和 截 距 不 变性 的 严格 限制 。 模 型 
对 应 的 Mplus 代码 见 De Bondt 和 Van Petegem (2015) 的 补充 材料 。 

多 组 模型 结果 如 下 : 过 度 激 动 问卷 在 五 个 维度 上 均 符 合 跨 组 截 距 不 变性 , 即 男女 性 两 组 
的 条 目 载荷 和 截 距 符合 跨 组 相同 , 模型 拟 合 良 好 (智力 维度 : PPp 值 为 0.540; 想象 维度 : PPp 
值 为 0.392; 情感 维度 : PPp 值 为 0.500; 感官 维度 : PPp 值 为 0.598; 精神 运动 维度 : PPp 
值 为 0.518)。 通 过 比较 组 间 因 子 均 值 的 差异 发 现 ， 女 性 的 因子 均值 在 情感 和 感官 维度 上 均 显 
著 高 于 男性 , 而 男性 的 因子 均值 则 在 精神 运动 维度 上 显著 高 于 女性 。 在 其 它 维 度 上 未 发 现 因 
子 均值 的 跨 性 别 差 异 。 


ANS 


3.4 贝 叶 斯 多 层 结构 方程 模型 


聚 类 抽样 (Cluster Sampling) 是 心理 学 研究 中 一 种 常用 的 抽样 方法 。 聚 类 抽样 指 按 某 种 标 
准将 总 体 划分 为 群 (组 )， 随 机 抽取 部 分 群 (组 )， 并 以 该 部 分 群 (组 ) 中 的 部 分 个 体 为 样本 。 而 这 
种 抽样 法 产生 的 数据 为 多 层 ( 撕 套 ) 数 据 结构 ， 如 从 不 同班 级 中 抽取 学 生 样 本 ， 其 中 ， 学 生 层 
次 幅 套 于 班级 层次 。 而 上 文 提 到 的 追踪 数据 也 可 作为 多 层 数据 中 的 一 种 , 其 中 不 同 的 测量 时 
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针对 该 类 数据 结构 ， WR RAEN TES BUR URE a, 而 且 会 违反 独立 
性 假设 , 因为 同一 群 (组 ) 内 的 变量 可 能 存在 相关 。 而 采用 多 层 模型 可 以 将 方差 分 解 到 各 层次 ， 
进而 探究 不 同 层 次 中 变量 间 的 关系 。 在 潜 变 量 建 模 中 ， 多 层 结构 方程 模型 (Multilevel 
Structural Equation Modeling, MSEM) 人 允许 模型 具有 来 自 组 间 和 组 内 的 不 同方 差 与 协 方差 ， 且 


可 以 处 理 超过 两 层 的 艇 套 数据 (Rabe-hesketh, Skrondal, & Pickles, 2004). 


但 由 于 多 层 数据 结构 和 建 模 的 复杂 性 , 传统 的 极 大 似 然 估计 容易 在 参数 估计 及 模型 收敛 
和 拟 合 上 出 现 问题 。 而 且 当 每 一 层 的 平均 样本 量 减 少时 , 参数 估计 准确 性 和 模型 的 收敛 率 会 


下 降 (Hox, Maas, & Brinkhuis, 2010; Li & Beretvas, 2013; Meuleman & Billiet 2009; Preacher, 


Zhang, & Zyphur, 2011)。 此 外 , 组 内 相关 系数 (Interclass Correlation Coefficient, ICC) 也 会 影响 
参数 估计 的 质量 ，ICC 指 组 间 变 异性 和 总 体 变异 性 的 比值 。 在 进行 MSEM 分 析 时 ， 如 果 同 
时 出 现 组 内 相关 系数 低 和 每 一 层 的 平均 样本 量 少 的 情况 , 容易 导致 协 方差 矩阵 不 正定 和 错误 
的 参数 估计 ， 如 ， 残 差 方 差 为 负数 (Depaoli & Clifton, 2015; Li & Beretvas, 2013). 

相 较 于 极 大 似 然 估计 ， 贝 叶 斯 估计 在 多 层 结构 方程 建 模 中 表现 更 好 : 首先 ， 贝 叶 斯 方法 
对 参数 估计 的 准确 性 更 高 (Asparouhov & Muthén, 2010; Baldwin & Fellingham, 2013); 其 次 ， 
对 于 参数 提供 一 定 的 先 验 信 息 也 可 以 避免 参数 估计 中 出 现 负 方 差 或 模型 无 法 收敛 的 问题 


(Chung, Rabe-Hesketh, Dorie, Gelman, & Liu, 2013); 此 外 , 贝 叶 斯 方法 更 易于 处 理 复杂 模型 ， 


如 三 层 且 含有 有 序 分 类 变量 数据 的 模型 (Asparouhov & Muthén, 2012); 最 后 ，MSEM 中 最 常 


见 的 问题 是 模型 第 二 层 的 样本 量 较 小 , 而 贝 叶 斯 方法 在 小 样本 情况 下 依旧 能 够 提供 准确 的 参 


数 估 计 (Asparouhov & Muthén, 2010; Baldwin & Fellingham, 2013). 


在 实证 研究 中 , Johnson, Schoot, Delmar 和 Crano (2015) 通 过 建立 贝 叶 斯 多 层 结构 方程 模 
型 ， 探 究 了 在 小 组 合作 的 两 个 人 际 交往 过 程 即 讨论 和 争论 之 间 的 相互 影响 和 促进 ;而 
Tamminen, Gaudreau, Mcewen 和 Crocker (2016) 为 了 避免 传统 方法 下 容易 出 现 模 型 无 法 收敛 
的 问题 ， 建 立 贝 叶 斯 MSEM 检验 影响 运动 员 享 受 程度 的 因素 ， 发 现 运动 员 个 人 情绪 管理 策 
略 、 和 群体 层面 的 队伍 气氛 和 同伴 关系 均 与 享受 程度 有 关 。 

为 了 方便 研究 者 应 用 该 方法 ,我 们 将 以 Prem, Scheel, Weigelt, Hoffmann 和 Korunka (2018) 
的 研究 为 例 , 展示 如 何 将 该 方法 应 用 到 实际 数据 分 析 。Prem 等 人 (2018) 采 用 了 贝 叶 斯 多 层 结 
构 方程 模型 来 探究 工作 特征 影响 工作 场所 中 的 拖延 的 中 介 机 制 。 该 研究 收集 了 110 名 雇员 的 
数据 ， 每 位 员工 需 完成 为 期 12 天 、 每 天 三 次 的 问卷 填写 。 

问卷 包括 四 个 变量 : 自 变量 即 工作 特征 , 包括 时 间 压 力 、 问 题解 决 和 计划 与 决策 三 个 维 
度 ; 两 个 中 介 变 量 即 认 知 评估 过 程 和 自我 调节 努力 , 其 中 认 知 评估 过 程 包括 挑战 评估 和 阻碍 
评估 两 个 维度 ; 因 变 量 即 工作 场所 拖延 ; 睡眠 质量 和 职业 自我 效能 感 作为 控制 变量 。 被 试 内 
和 被 试 间 两 个 水 平 的 假设 一 致 , 均 如 图 3 所 示 : 假设 1 是 工作 特征 通过 增加 挑战 评估 来 减少 
自我 调节 努力 从 而 减少 影响 工作 场所 拖延 ; 假设 2 是 工作 特征 通过 增加 阻碍 评估 来 增加 自我 
调节 努力 从 而 增加 工作 场所 拖延 。 

由 于 数据 为 嵌 套 结构 ,因此 研究 者 采用 多 层 结构 方程 模型 进行 分 析 , 将 变量 方差 分 解 至 


被 试 间 部 分 和 被 试 内 部 分 。 相 较 于 其 他 多 层 中 介 分 析 , 结构 方程 建 模 产生 的 偏差 更 小 。 此 外 ， 
大 多 数 情 况 下 被 试 内 水 平 的 间接 效应 呈现 偏 态 分 布 , 因此 采用 贝 叶 斯 估计 还 可 以 更 好 的 检验 


间接 效应 (Preacher, Zyphur, & Zhang, 2010)。 

研究 者 采用 Mplus 8.0 进行 贝 叶 斯 多 层 结构 方程 模型 分 析 , 结果 表明 模型 拟 合 良好 , PPp 
值 为 0.460。 参 数 估计 结果 表明 : 被 试 内 水 平 中 ， 工 作 特 征 的 三 个 维度 均 对 工作 场所 拖延 有 
负 向 的 序列 间接 效应 , 即 通过 增加 挑战 评估 来 减少 自我 调节 努力 , 从 而 影响 工作 场所 拖延 (时 
间 压 力 : 间接 效应 大 小 为 -0.004，95%CI[-0.009，-0.000]; 问题 解决 : 间接 效应 大 小 为 -0.011， 
95%CI[-0.023, -0.002]; 计划 与 决策 : 间接 效应 大 小 为 -0.004, 95%CI[-0.009, -0.000])。 此 外 ， 


时 间 压 力 对 工作 场所 拖延 还 有 正 向 的 序列 间接 效应 , 即 通过 增加 阻碍 评估 来 增加 自我 调节 努 


力 从 而 影响 工作 场所 拖延 (间接 效应 大 小 为 0.006, 95%CI[0.002, 0.012])。 而 在 被 试 间 水 平 中， 


上 述 的 序列 间接 效应 并 不 成 立 。 因此, 该 研究 揭示 了 工作 特征 与 工作 场所 拖延 存在 密切 关系 ， 


且 主 要 通过 被 试 内 水 平 上 的 认 知 评估 和 自我 调节 努力 来 实现 的 。 


工作 特征 认 知 评估 过 程 自我 调节 努力 工作 拖延 


— |) a 
[saa |p 
— a 
计划 与 决策 


图 3 研究 概念 模型 图 


4 ”模型 评价 与 拟 合 指 标 


传统 的 模型 拟 合 指标 并 不 适用 于 贝 叶 斯 估计 , 因此 在 使 用 贝 叶 斯 方法 时 , 需要 用 到 以 下 


指标 对 模型 进行 有 效 的 评估 。 


4.1 后 验 预 测 p 值 


模型 和 数据 的 拟 合 程度 可 以 通过 后 验 预 测 检 验 来 进行 评估 (Gelman， Meng, Stern, & 
Rubin，1996)。 后 验 预测 检验 比较 了 实际 数据 与 假设 模型 产生 的 数据 之 间 的 差异 ， 可 以 用 于 
评估 模型 和 实际 数据 的 拟 合 程度 。 通 过 后 验 预测 检验 可 以 得 到 后 验 预 测 p 值 (Posterior 


Predictive p-value, PPp 值 )。PPp 值 与 假设 检验 中 的 p EELDE, CE MCMC 算法 多 


长 


次 迭代 中 ， 依 据 理论 模 型 生成 的 统计 检验 量 大 于 样本 数据 的 统计 检验 量 的 比例 。 因 此 ，PPp 
值 在 0.5 左右 ， 即 接近 随机 概率 1/2 时 ， 表 示 模 型 拟 合 得 很 好 。 此 外 ， 后 验 预 测 检验 还 能 够 


给 出 样本 数据 与 模 
数 ， 且 0 落 在 区 间 


4.2 贝 叶 斯 因子 


型 生成 数据 之 间 统 计 检 验 量 差异 的 95% 置 信 区 间 ， 当 该 区 间 的 下 限 为 负 


中 心 时 ， 表 示 模 型 拟 合 得 很 好 (Muthkn & Asparouhov, 2012). 


贝 叶 斯 因子 (Bayes FactonD 是 用 于 模型 比较 的 重要 统计 量 ， 它 可 用 于 比较 非 花 套 模 型 。 它 
可 以 比较 在 已 有 数据 集 下 , 支持 Mo 与 Mi 两 个 模型 的 概率 。 这 样 的 比较 不 依赖 于 假设 检验 ， 
即使 样本 量 很 大 也 不 会 倾向 于 支持 备 择 假 设 Mi( 李 锡 钦 ，2011)。Kass 和 Raftery (1995) 提 出 


了 解释 贝 叶 斯 因子 


的 准则 。 如 ， 贝 叶 斯 因子 介 于 1 到 3 之 间 表 示 该 数据 对 两 个 模型 的 支持 程 


度 差不多 ， 此 时 在 模型 选择 时 还 需要 考虑 “简约 性 ”原则 ， 或 结合 其 他 指标 进行 判断 。 


4.3 贝 叶 斯 信息 准则 


贝 叶 斯 信息 准 


则 (Bayesian Information Criterion, BIC; Schwarz，1978) 与 赤 池 信息 准则 


(Akaike Information Criterion, AIC) 一 样 ， 在 评价 模型 拟 合 的 同时 考虑 了 模型 复杂 度 ， 对 复杂 


的 模型 进行 惩罚 ， 


拟 合 更 好 。 但 BIC 是 相对 拟 合 指 标 ， 不 能 评价 两 个 模型 拟 合 差异 的 显著 性 。 


4.4 偏差 信息 准则 


且 BIC 比 起 AIC 对 复杂 模型 的 惩罚 作用 更 强 。BIC 的 值 越 小 ， 说 明 该 模型 


与 BIC 一 样 ，1 


(Spiegelhalter, Best, 


ii Ze (a JS YE UU] (Deviance Information Criterion, DIC) 也 被 用 来 比较 竞争 模型 


Carlin, & Van der Linde, 2002)。DIC 越 小 表示 模型 拟 合 越 好 。 但 是 由 于 DIC 


更 加 符合 贝 叶 斯 偏差 (Bayesian Deviance) 的 概念 (Kaplan & Depaoli, 2012)， 其 应 用 更 为 普遍 。 


相 比 之 下 ，BIC 更 多 被 用 于 频率 学 方法 中 的 模型 比较 。 


5 软件 介绍 


尽管 采用 贝 叶 


斯 方法 估计 SEM 有 着 诸多 优点 ， 但 其 应 用 一 直 较 为 滞后 ， 一 个 重要 的 原 


因 就 是 它 看 起 来 要 求 研 究 者 具备 很 好 的 贝 叶 斯 统计 基础 (MuthEn & Asparouhov, 2012)。 但 实 


际 上 目前 能 够 采用 


于 学 习 使 用 。 其 中 


贝 叶 斯 方法 分 析 SEM 的 软件 不 仅 可 以 满足 大 多 数 研究 者 的 需求 ， 而 且 易 


常用 软件 主要 有 以 下 几 个 ， 除 Mplus 外 其 余 软 件 均 为 免费 开源 软件 : 


(1) 目前 最 为 流行 的 潜 变 量 分 析 软 件 Mplus (Muthén & Muthén, 1998~2017), 其 编程 语言 


简单 易学 , 由 于 其 默认 设 定 较 多 , 在 使 用 中 研究 者 一 般 只 需要 将 估计 方法 从 极 大 似 然 估计 换 
为 贝 叶 斯 佑 计 ， 并 提供 相应 先 验 信息 即 可 。 如 果 没 有 提供 先 验 信息 Mplus 通常 会 默认 提供 
无 信息 先 验 分 布 。 不 过 Mplus 8.0 还 不 能 估计 一 些 特殊 模型 ， 如 潜 调 节 模 型 ( 即 调节 变量 为 潜 
变量 的 模型 ) 等 。 当 研究 者 将 估计 方法 设 定 为 贝 叶 斯 估计 时 ，Mplus 8.0 能 提供 PPp 值 、DIC 


和 BIC 值 作为 模型 拟 合 与 评价 的 指标 ; 


(2) WinBUGS (Windows version of Bayesian Inference Using Gibbs Sampling; Lunn et al., 
2000) 是 专门 用 于 贝 叶 斯 统计 推断 的 软件 包 。 相 比 于 Mplus， 其 对 研究 者 的 贝 叶 斯 统计 基础 
有 更 高 的 要 求 。Mplus 中 一 些 默 认 的 估计 设 定 在 WinBUGS 中 需要 研究 者 自己 设 定 ， 但 是 它 
功能 强大 ， 可 以 灵活 地 对 复杂 模型 进行 估计 。WinBUGS 能 提供 DIC 值 作为 模型 拟 合 与 评价 
的 指标 ; 

(3) 软件 Stan (Stan Development Team, 2014) 可 以 对 复杂 的 BSEM( 如 潜 调 节 模 型 、 多 层 
模型 等 ) 进 行 估计 ， 且 可 以 与 最 流行 的 数据 分 析 语 言 (如 R, Matlab, Python 等 ) 接 口 ; 


(4) R 软件 中 的 blavaan (Merkle, & Rosseel, 2015) 尚 不 能 估计 一 些 特殊 的 模型 ,如 潜 调 


节 模 型 、 含 有 序 分 类 变量 的 模型 等 。 不 过 用 户 可 以 根据 需要 导出 JAGS (Just Another Gibbs 
Sampler; Plummer, 2005) 代 码 来 估计 复杂 、 特 殊 的 模型 。 此 外 ， 通 过 mplus2lavaan0 函 数 还 可 
以 将 Mplus 软件 与 blavaan 接口 。blavaan 能 提供 PPp 值 以 及 DIC 值 作为 模型 拟 合 与 评价 的 
指标 。 

在 本 文 所 介绍 的 应 用 研究 中 ， 大 部 分 研究 都 采用 了 Mplus 软件 进行 贝 叶 斯 结构 方程 建 


模 (例如 ,Crenshaw et al., 2016; Falkenstrom et al., 2015; Golay et al., 2013; Zeman et al., 2017)， 


但 也 有 两 篇 研究 采用 R 和 JAGS 软件 进行 建 模 (Praetorius et al., 2017; Winans-Mitrik et al., 


2014)。 研 究 者 可 以 根据 自己 的 建 模 需要 选择 相应 的 分 析 软 件 。 


6 讨论 


在 结构 方程 建 模 中 贝 叶 斯 方法 有 着 无 可 替代 的 优势 ,在 模型 识别 和 拟 合 ， 参 数 估 计 ， 处 
里 复杂 模型 和 小 样本 情况 等 方面 贝 叶 斯 方法 都 有 着 更 好 的 表现 。 该 方法 在 近 几 年 的 发 展 也 颇 
为 迅 独 (Van de Schoot et al., 2017)。 基 于 贝 叶 斯 方法 结合 先 验 信息 ， 不 依赖 方差 协 方差 矩阵 
等 特性 ， 衍 生出 了 新 的 建 模 思 路 ， 如 放宽 对 传统 测量 模型 的 严格 限制 (Lu et al., 2016; Muthén 
& Asparouhov, 2012; Pan et al., 2017)， 多 组 模型 中 允许 组 间 存 在 较 小 的 差异 (Muthén & 


Asparouhov, 2013) 等 ， 而 这 些 建 模 思 路 在 传统 方法 中 都 难以 实现 。 


ie 


chinaXiv:201812.00865v1 


但 是 截至 2018 年 6 月 ， 在 中 国 知 网 数据 库 检索 发 现 ， 在 国内 心理 学 期 刊 中 尚未 有 贝 叶 斯 


结构 方程 模型 的 应 用 研究 。 和 希望 本 文 可 以 为 国内 心理 
在 应 对 小 样本 、 传 统 方法 下 模型 无 法 识别 等 情况 时 可 以 游 胃 有 余 。 此 外 ， 频率 学 派 方法 的 主 


导 地 位 使 得 心理 学 研究 者 对 贝 叶 斯 方法 的 了 解 不 足 。 
斯 的 了 解 是 相当 有 限 的 , 这 使 得 应 用 研究 者 在 面 对 这 种 新 的 研究 工具 时 可 能 会 担心 无 法 很 好 


地 掌握 这 种 方法 。 但 实际 上 对 于 有 结构 方程 


EE 模 基 而 


学 研究 者 在 实证 研究 中 提供 新 的 思路 ， 


王 孟 成 等 人 (2017) 也 指出 学 术 界 对 贝 叶 


的 研究 者 来 说 ， 采 用 Mplus 软 件 进行 贝 


叶 斯 估计 是 非常 易于 掌握 的 ， 希 望 本 文 的 介绍 可 以 打破 这 种 刻板 印象 。 
尽管 使 用 贝 叶 斯 方法 有 着 诸多 优势 , 但 先 验 的 潜在 影响 、 贝 叶 斯 特征 与 结果 的 错误 解读 ， 
以 及 不 正确 的 贝 叶 斯 结果 报告 都 可 能 和 带 来 错误 。 针 对 以 上 几 点 潜在 的 危险 , Depaoli 和 van de 


Schoot (2015) 提 出 了 WAMBS 清 单 以 避免 贝 叶 斯 方法 的 误 用 ， 研 究 者 可 以 参考 这 个 清单 以 避 
免 结果 报告 的 错误 。 而 随 着 采用 贝 叶 斯 方法 估计 结构 方程 模型 的 应 用 和 方法 研究 越 来 越 多 


(Van de Schoot et al., 2017)， 相 信 其 建 模 方法 、 使 用 技巧 和 报告 规范 会 越 来 越 完 善 。 


NI 
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Bayesian structural equation modeling and its current researches 
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‘a Department of Psychology, Sun Yat-sen University, Guangzhou 510006, China) 
Abstract: Structural equation modeling (SEM) has been widely used in psychological researches 
to investigate the casual relationship among latent variables. Model estimation can be conducted 
under both the frequentist framework (e.g., maximum-likelihood approach) and the Bayesian 
framework. In recent years, with the prevalence of Bayesian statistics and its advantages in 
dealing with small samples, missing data and complex models in SEM, Bayesian structural 
equation modeling (BSEM) has developed rapidly. However, in China its application in the field 
of psychology is still insufficient. Therefore, this paper mainly focuses on presenting this new 
research method to applied researchers. We explain the theoretical and methodological basis of 
BSEM, as well as its advantages and disadvantages compared with the traditional frequentist 
approach. We also introduce several commonly used BSEM models and their applications. 
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